High-energy excitonic effects in graphite and graphene 
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We present ah initio many-body calculations of the optical absorption in bulk graphite, graphene 
and bilayer of graphene. Electron-hole interaction is included solving the Bethe-Salpeter equation 
on top of a GW quasiparticle electronic structure. For all three systems, we observe strong excitonic 
effects at high energy, well beyond the continuum of tt — s- tt* transitions. In graphite, these affect the 
onset of 0" — ^ (J* transitions. In graphene, we predict an excitonic resonance at 8.3 eV arising from 
a background continuum of dipole forbidden transitions. In the graphene bilayer, the resonance is 
shifted to 9.6 eV. Our results for graphite are in good agreement with experiments. 

PACS numbers: 71.35.-y, 71.35. Cc, 78.20.Bh, 78.40.-q 
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Graphene, a recently discovered 2D hexagonal crys- 
tal carbon sheet [ij, has attracted much interest due to 
its exotic electronic properties. The comparison with 
ordinary 3D graphite, the ABA stacking of graphene 
layers with weak inter-layer interactions, gives further 
insights, as screening effects and effective 2D confine- 
ment are modified there. Both systems share a peculiar 
semi-metallic character together with a strong electronic 
anisotropy, giving rise to optical properties of particular 
interest, especially in view of technological applications in 
opto-electronics More generally, also in astrophysics, 
accurate determination of optical absorption of carbon 
structures is fundamental 

The optical absorption of CTaphite was experimentally 
determined by reflectance [3Jl| as well as by energy-loss 
via Kramers-Kronig , showing important anisotropy 
effects between measurements using polarizations paral- 
lel or perpendicular to the crystallographic c axis. Op- 
tical properties of graphene are under current experi- 
mental investigation |1CH12| . A first theoretical JDOS 
(joint density of states) calculation based on an in- 
dependent particle picture and usual selection rules, has 
shown, that the optical absorption spectrum of graphite 
can be divided in two regions: the visible range from 
eV up to 5 eV originates from transitions among tt bands, 
whereas the region beyond 10 eV is made of a band tran- 
sitions. Ab initio calculations beyond JDOS and in the 
RPA approximation have been done on graphite 14| and 
graphene, also including local- field effects |15j . 

In this work, we extend previous theoretical results 
[3, [l^ beyond RPA by including electron-electron ( e- 
e) and electron-hole (e-/i) interactions. Our calcula- 
tions are based on the many-body ab initio GW and 
Bethe-Salpeter equation (BSE) approach [l^. We study 
bulk graphite, free-standing graphene, and a bilayer of 
graphene, considering both polarizations. With respect 
to Ref. [l6| which is restricted to the low-energy range, 
our work extends also to high energies. 



We observe strong excitonic effects at unusual high en- 
ergies, well beyond the continuum of tt-tt* transitions. In 
graphite, excitonic effects strengthen the low energy part 
of the a-a* structure and reshape it. Both, in the perpen- 
dicular and in the parallel polarization, our spectra are in 
very good agreement with experiment, in particular with 
Ref. [8|. In graphene, we predict an intense peak at 8.3 eV 
that we interpret as an excitonic resonance arising from 
a background single-particle continuum of dipole forbid- 
den transitions. In the bilayer the resonance remains, 
but it is shifted to 9.6 eV due to reduced confinement 
and increased screening. Measuring optical spectra, this 
feature could thus be used as a fingerprint to discriminate 
between graphene and multi-layer graphene. 

Our GW and BSE calculations [l7| are based on 
ground state calculations using density-functional theory 
in the local-density approximation (DFT-LDA). Starting 
from the Kohn-Sham DFT electronic structure, we cal- 
culate an ab initio GW quasiparticle electronic structure 
that takes into account e-e many-body interactions. In 
order to include e-h interactions in the response func- 
tions, we solve the Bethe-Salpeter equation for the two- 
particle correlation function L, 



L = GG + GGEL, 



(1) 



where G is the GW Green's function, and S is the BSE 
kernel. We have used the approximation S — —iv + iW , 
where v is the Coulomb and W the screened inter- 
action. As described in Ref. [13], we define H'^^'^ = 
(e^w _ gGW) _,_ 2^ and remap the BSE, Eq. (P), into a 
2-particle Schrodinger equation 
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where 
ing the 



J^CXC 



represent the excitation energies includ- 
e-h interaction effects, and '^'^^ the exci- 
tonic wavefunctions. We diagonalize Eq. ([2]) working 



with a basis of Kohn-Sham bilinears, '^"^'^ 



2 




graphite 



out-of-plane 





EXP TB 




EXPZ 




EXPK 




EXP V 




KS-RPA 




KS-RPA NLF 




GW-RPA 




BSE 




(i)[cVl 



(B[cVl 



FIG. 1: (color online) Optical absorption spectra of bulk 
graphite for E _L c. Solid black line: BSE; blue dashed 
line: GW-RPA; green dot-dashed line: KS-RPA. All theo- 
retical curves convoluted by a relative Gaussian broadening 
of o" = 0.075 a;. Experiments: TP, cyan triangles down [3|; 
TB, orange diamonds: 0; Z, indigo squares: [g|; K, magenta 
triangles up: 0]; V, red circles: [Ij. 



J2kvc'^\'^'l^vk i'''h)4'cki'''e) where V (c) runs on valence 
(conduction) bands, and k lies in the 1st Brillouin zone. 
The macroscopic dielectric function is then given by 
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The ground-state DFT-LDA and the GW corrections 
have been calculated using the ABINIT code. In the case 
of the bilayer and of graphene, the layers are isolated by 
38 Bohr of vacuum, distance large enough to avoid spu- 
rious interactions between replicas. We used Martins- 
Trouiller pseudopotentials with s and p electrons in the 
valence. The BSE calculation was carried out by the EXC 
code. For graphite, the Brillouin zone was sampled with 
a (8 8 5) Monkhorst-Pack k-point grid shifted by (0.01, 
0.02, 0.03), whereas a (16 16 1) shifted grid was used 
for both graphene and bilayer. Wavefunctions have been 
represented using 967 plane waves (39 Ry) for graphite 
and 1367 (24 Ry) for the bilayer and graphene, while the 
dimension of the kernel S was 53 (graphite) and 287 (bi- 
layer and graphene) plane waves (8 Ry). We included 
all cr and tt occupied bands, i.e. 8 in graphite and the 
bilayer, 4 in graphene. And 13 empty bands in graphite, 
25 in the bilayer and 21 in graphene. We also used the DP 
code to obtain random-phase approximation (RPA) spec- 
tra. Local- field effects have been taken into account, both 
in RPA, and in BSE calculations. All theoretical curves 
have been convoluted by a Lorentzian with a quadratic 
width rj = 0.003 (Fermi liquid-like lifetime behavior) 
to introduce the intrinsic quantum broadening, together 
with a Gaussian broadening adapted for the presumed 
experimental energy resolution. 



FIG. 2: (color online) Optical absorption spectra of graphite 
for E II c. Same notation as in Fig. [1] but with a relative 
Gaussian broadening of a = 0.025 oj. 



Graphite, in-plane polarization: The calculated imagi- 
nary part of the macroscopic dielectric function, directly 
associated with the absorption spectra, is shown in Fig. [T] 
for graphite in the E _L c polarization. The Kohn-Sham 
RPA (KS-RPA) spectra reproduces previous theoretical 
results [llj calculated using a denser grid of k-points. 
Local-field (LF) effects are expected to be small since 
the in-plane electronic density of graphene is almost ho- 
mogeneous [15i] . Indeed, spectra with and without LF 
coincide for this polarization. The position of the lowest 
energy 4.2 eV peak of the KS-RPA spectra is exactly the 
same as in all the considered experiments. On the other 
hand, the position of the highest energy structure due to 
a ^ a* \s compatible with only some of the experiments. 

The inclusion of e-e GW effects globaly shifts the spec- 
trum to higher energies. The shift however is not rigid. 
GW corrections are less effective on bands closer to the 
Fermi energy around K, but they increase with energy, 
and correct the band-gap underestimation of the DFT- 
LDA electronic structure, for example at F and L, re- 
sulting in an increased agreement with photoemission ex- 
periments. However, concerning optical absorption, the 
agreement with the experiment (see Fig. [1]) gets worse 
when passing from KS-RPA to GW. For example, the 
KS-RPA main structure at 14 eV has been shifted to 15.3 
eV by GW corrections, at least 1.3 eV off from the ex- 
perimental position, and, similar, also the lowest energy 
structure around 4 eV, shifted by 0.7 eV. 

The inclusion of e-h interaction effects via BSE seems 
to compensate e-e effects, restoring a good agreement 
with experiment, similar to the KS-RPA result. The low 
energy peak position agrees with all five considered ex- 
periments, although the height and the shape differ some- 
how. The position of the high energy peak recovers a po- 
sition at smaller energy than the KS-RPA peak. Its pro- 
file is slightly reshaped, with a strengthening of the peak. 
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FIG. 3: (color online) Optical absorption spectrum of 
graphene for E _L c. Inset for E || c. Same notation as in 
Fig. [T]witli an absolute Gaussian broadening of 0.2 eV. 



FIG. 4: (color online) Optical absorption spectrum of bilayer 
of graphene for E _L c. Same notation as in Fig. [T] with an 
absolute Gaussian broadening of 0.2 eV. 



particulary on the low-energy side. With respect to KS 
or GW-RPA, the peak position appears now to be at 12.8 
eV. A fine analysis reveals a main excitonic energy at 12.6 
eV, together with two other main excitation energies at 
13.3 and at 13.7 eV. The latter should conjure an asym- 
metric, slower drop on the right shoulder of the peak. 
Unfortunatly, there is some disagreement between exper- 
imental spectra in this high energy region measured by 
different techniques mostly energy-loss and reflec- 

tivity measurements, and the position of the high-energy 
peak varies between 12.6 and 14.3 eV. The position of the 
peak of our BSE calculation is more in agreement with 
the results of Venghaus Q and also Tosatti and Bassani 
0, both at 12.6 eV and both measured by energy-loss, 
which do not suffer from surface effects that tamper the 
bulk result for reflection experiments. The shoulder ob- 
served at '--^ 14 eV could correspond to our excitation at 
13.7 eV, although the imposed broadening (linear Gaus- 
sian with a = 0.075 a;) makes it less evident in the figure. 
The agreement of our result with the optical experiment 
of Taft and Philipp Q is still acceptable, but worsens 
when comparing with Zeppenfeld [dl (peak at 14 eV) and 
in particular with Klucker et al. [5| (peak at 14.3 eV), 
measured via X-ray absorption. Comparing BSE with 
GW, it is remarkable that excitonic effects are stronger 
in the high-energy range of the spectrum rather than in 
the lower one, in contrast to ordinary semiconductors or 
insulators where excitonic effects mainly affect the low- 
est energy part of optical spectra, or to metals, where 
excitonic effects are in general negligible. 

Graphite, out-of-plane polarization: Previous conclu- 
sions are further confirmed by the results for the E || c 
polarization (Fig. [5]). As already found in Ref. [l^, LF 
effects are stronger in the out-of-plane case, as it is evi- 
dent comparing Kohn-Sham RPA with and without LF. 
Both curves are in agreement with previous theoretical 



results IJ, [ll]. The KS-RPA with LF curve IS again m 
fair agreement with the experiment. The position of the 
3.7 eV low energy peak is quite accurate, the main peak 
at 11.5 eV not far from the experiment, the last peak at 
16.3 eV and a shoulder at ^ 13.2 eV in correspondence 
with experimental features. Again, the inclusion of e-e 
GW interaction effects causes a blue shift of up to 1 eV 
and worsens the agreement with the experiment. The im- 
provement appears only when including also e-h interac- 
tion effects (BSE) . The peak at low energy is now between 
the KS-RPA and the GW-RPA peaks, at exactly 4.2 eV, 
in agreement with experiment, especially with Venghaus 
iSj. The main peak instead red-shifts beyond the KS- 
RPA by 0.6 eV and is placed at 10.9 eV, again in good 
agreement with Venghaus (10.8 eV) and at lower energy 
with respect to Zeppenfeld [63 and Tosatti and Bassani 
y] (both at 11.3 eV). With respect to KS-RPA and GW, 
e-h effects also enhance the peak, and improve the agree- 
ment with the experiment on the magnitude. The higher 
energy peak gets restored at the KS-RPA position with 
a slight improvement on the position, but a worsening 
in intensity. The ~ 13 eV shoulder, still evident in KS- 
RPA, has been washed out by the strengthening of the 
neighbor main peak and is no longer visible. Our ac- 
tual BSE calculation supports the scenario supposed for 
many-body effects by Ref. [l5l |. 

With respect to experiment, in the case of out-of-plane 
polarization, the situation is similar to the in-plane po- 
larization. The BSE result is in particular in good agree- 
ment with Venghaus Q, but also with Tosatti and Bas- 
sani Hand Zeppenfeld Q, but disagrees with Klucker 
et al. [5j- However, apart from the peak positions, the 
agreement with the experiment can be improved concern- 
ing peak intensities. In particular the low-energy part of 
the spectrum, due to interlayer interactions, still seems 
to be underestimated. There is also a slight underestima- 
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tion on the main peak, although for this polarization we 
used less broadening (cr = 0.025 a;) than in the in-plane 
polarization. 

Graphene: In graphene (Fig. [3]), due to its effective 2D 
dimensionality, LF effects become the most important ef- 
fects for the out-of-plane polarization. As already found 
in Ref. , for E || c, LF effects produce a complete sup- 
pression of spectra due to a depolarization effect (com- 
pare KS-RPA with and without LF curves in the inset of 
Fig.E]). The 2D selection rules forbid tt -J> tt* transi- 
tions (<10 eV), allowing a ^ a* ones instead. Neverthe- 
less, the strong anisotropy increases the LF effects sup- 
pressing the latter transitions as well. On the other hand, 
LF effects are negligible in-plane (Fig. ^ . In graphene as 
in graphite, GW e-e interaction effects produce a ~ 1 eV 
blue shift of the KS spectrum. Both the tt — )■ tt* peak at 
4 eV and the cr — >■ a* peak at 14 eV are shifted around 
+1 cV. In graphene, however, the inclusion of e-h interac- 
tion turns out to cause a much more spectacular effect. 
The spectrum is completely reshaped. Surprisingly, at 
4^5 eV, as well as at 14 ~ 15 eV, there is not anymore 
any oscillator strength. The RPA peaks disappear com- 
pletely. Most of the spectral weight is pushed to very 
low energies, < 3 eV. However, unexpectedly, a strong 
peak rises up at 8.3 eV in the region where single-particle 
oscillator strengths were vanishing. This excitonic reso- 
nance arises from a background single-particle continuum 
of dipole forbidden transitions. An analysis over the ex- 
citonic oscillator strengths |5'^'^'^(wfc -I- q\e~^'^^\ck)\ rela- 
tive to the most intense excitonic eigenvalue indicates 
that the strong excitonic resonance is conjured by an ex- 
tended mixing of cr —> cr* transitions over a large en- 
ergy range, from 12.5 to 17 eV. An excitonic effect at so 
large distance in energy (8.3 vs 12.5-17 eV) is caused by 
the increased e-h interaction due to reduced screening in 
graphene. Since we do not have any absorption spectrum 
experiments in graphene for comparison, we maintained 
the broadening at the minimum level to highlight these 
effects. Depending on the experimental energy resolu- 
tion, this peak could appear more broadened, but does 
not affect the peak position at < 8.5 eV, and a slight 
right asymmetry. We note that inclusion of band 2 in 
the BSE is the only crucial ingredient to start to conjure 
the high energy excitonic resonance, whereas the 4 5 
eV |16| energy region is very sensitive to the total num- 
ber of bands, the dimension of the BSE kernel, and the 
vacuum distance between graphene layers. 

Bilayer: Similar to graphene, an excitonic resonance 
also occurs for the graphene bilayer (Fig. 2]), however, 
at larger energy, 9.6 eV. The difference is due to the 



increased screening with respect to graphene, which re- 
duces e-h interaction strength. With respect to graphene, 
RPA spectra appear different at low energy due to in- 
terlayer effects. However, also here, e-h interaction ef- 
fects sweep the region at ~ 5 eV, and push the oscillator 
strength towards lower energies < 3 cV, or towards the 
excitonic resonance. 

In conclusions, we have presented ah initio GW-BSE 
calculations of the optical absorption spectra in graphite, 
graphene, and the bilayer, taking into account e-e, e-h, 
and LF effects. For the three systems, we have found 
strong excitonic effects in a high-energy range > 8 eV, 
well beyond the range of tt — > tt* transitions. In graphene 
and bilayer, these give rise to an excitonic resonance on 
a background continuum of dipole forbidden transitions. 
The excitonic resonance has a main a ^ a* character. 
Our BSE spectra on graphite are in good agreement with 
experiment. 
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